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ABSTRACT 

The evolution ot the contact interface and 
the propagation of compression waves inside a 
single wave rotor passage gradually opening to 
and traversing an inlet port is studied numerically 
using an inviscid formulation of the governing 
equations. Insights into the response of the 
interface and kinematics of the flow field to 
various opening times are given. Since the 
opening time is inversely proportional to the 
rotational speed of the rotor, the effects of 
passage rotation such as centripetal and Coriolis 
accelerations are intrinsically coupled to the 
gradual opening process. Certain three- 
dimensional features associated with the gradual 
opening process as a result of centripetal and 
Coriolis accelerations are illustrated. For the 
range of opening times or rotational speeds 
considered, a portion of the interface behaves 
like a vortex sheet that can degenerate into a 
complex interfacial structure. The vortices 
produced along the interface can serve as a 
stirring mechanism to promote local mixing. 
Coriolis and centripetal accelerations can 
introduce three-dimensional effects such as 
interfacial distortions in meridional planes and 
spanwise migration of fluid elements. 
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speed of sound 

dimensbniess ratio of time scales 
unit vector 

total energy per unit mass 
rothalpy 
passage length 
wheel Mach number 
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pressure 
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hub radius 
tip radius 
temperature 
time 

reference time scale 
relative velocity vector 
ratio of specific heats 
density 

azimuthal coordinate in a 
stationary reference frame 
azimuthal coordinate on the 
rotating passage 
passage sector angle 
rotational speed 


Subscripts 

L 

0 

oo 

Q 


longitudinal direction 
passage opening 
uniform reference state 
axis of rotation 


INTRODUCTION 

There is strong impetus for increasing the 
maximum cycle temperature of gas turbine 
engines in order to achieve higher cycle 
efficiencies. A wave rotor augmented cycle 
provides the potential for this increased 
performance within the limits of current materials 
technology [1]. Conceptually, the wave rotor is 
based on the idea of energy exchange using 
unsteady waves. A wave rotor cycle can play the 
role of a mechanically coupled compressor and 
turbine, delivering work to compress one gas 
stream while extracting work from another 
expanding gas stream. Since both expansion and 
compression are performed cyclicly in the rotor, 
the allowable material thermal limits are usually 
not exceeded. 

A wave rotor (Fig. 1) typically consists of 
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a number of cylindrical passages mounted along 
the periphery of a rotating drum. Each of these 
passages is cyclicly exposed to different end 
conditions that generate unsteady waves for 
energy transfer. As a wave rotor passage rotates 
into a stationary inlet-port, the passage is 
gradually exposed to the port flow. This process 
is called gradual opening and the period of 
adjustment wherein the leading and trailing wails 
of a rotor passage move out of a solid endwall 
allowing the passage to be fully exposed to the 
inlet-port conditions is herein referred to as the 
opening time. 

Eidelman [2] presented a numerical 
simulation of gradual opening in wave rotor 
passages showing that the dynamics of the 
passage opening significantly impact the flow 
pattern. The gradual opening process results in 
a significant distortion (i.e., obliquity) of the fluid 
interface in the passage-to-passage plane and 
also leads to a finite time for the formation of the 
shock. These events can adversely affect the 
performance of a wave rotor. Since Eidelman’s 
model was two-dimensional it could not include 
the effects of Coriolis and centripetal 
accelerations brought on by the rotation of the 
passage. 

For a given rotor passage geometry, the 
gradual opening time is inversely proportional to 
the rotational speed of the passage. Hence, the 
problem of gradual opening is intrinsically 
coupled to that of rotation. Larosiliere and Mawid 
[3] have shown that centrifugal and Coriolis 
accelerations brought on by the rotation of a 
channel can significantly impact unsteady 
gasdynamic wave processes. Therefore, the 
response of the flow in a wave rotor passage to 
the gradual opening process may be altered by 
Coriolis and centripetal accelerations resulting 
from the rotation of the passage. 

In the present study, the impact of 
gradual opening on the flow pattern in a wave 
rotor passage is re-examined in order to illustrate 
certain features due to Coriolis and centripetal 
accelerations that have not heretofore received 
much attention. A three-dimensional Euler model 
of the gradual opening problem is formulated in 
a rotating reference frame for a single wave rotor 
passage. The equations are numerically solved in 


order to track the gradual opening transient and 
the subsequent flow development as the passage 
traverses the inlet port. Variations in the 
response of the flow field in a passage of fixed 
geometry are examined for a wide range of 
rotational speeds. A discussion of the results 
highlighting implications for wave rotor design is 
provided. 


FORMULATION 

The gradual opening problem involves a 
coupling between the inlet-port flow and the flow 
in the wave rotor passages. Flow through the 
port is for the most part steady [1] , but the flow 
in the wave rotor passages is unsteady and 
periodic in time. This periodicity is a consequence 
of the cyclic nature of the wave rotor processes, 
interaction of the port flow with that in the 
passages during the gradual opening time may 
generate a local unsteady flow in the inlet port. 
This is probably a small scale unsteadiness since 
the width of the inlet port is usually much larger 
than that of a single passage. 

There are currently two different 
approaches for treating the gradual opening 
problem in a wave rotor. First, the more general 
approach adopted by Welch [4] involves the 
determination of a time-periodic flow field whose 
physical domain includes the inlet port and all of 
the passages spanning that port. The second 
approach called the Signal-Storage-and- 
Reconstmction (SSR) technique, introduced by 
Erdos et al. [5], determines the entire flow field by 
considering only the port and a single passage. 
The time periodicity condition is explicitly used to 
construct the flow in the other passages. Both 
techniques have been applied in rotor/stator 
interaction studies and are known to be 
computationally intensive in three-dimensions. 

In the present study, the aim is not to 
provide a detailed description of the coupling 
between the inlet-port flow and the flow in the 
rotor passage, but rather to characterize the 
response of the flow field in a rotating passage to 
the gradual application of the stationary inlet-port 
conditions. Thus, a single wave rotor passage is 
tracked as it passes stationary endwalls and the 
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inlet port. The flow in the inlet port is not 

calculated; however, its effect is accounted for =0 (1a) 

through boundary conditions at the inflow plane. 


Governing Equations 

The geometric configuration chosen to 
simulate the gradual opening problem is 
illustrated in Figure 2. It consists of a single 
cylindrical wave rotor passage of length L, and 
sector angle «I>. At t = 0 the passage is fully 
enclosed. Subsequent to this time, the solid-wall 
boundary conditions at the left side (x=0) of the 
passage are gradually replaced with those of the 
inflow-plane as the passage rotates out of the 
endwall. The right side is always closed. Viscous 
effects and heat conduction are neglected. 


B ^ + 2R^p§axW - 

(1b) 

= -^P 


B ^ + ^-(p =0 (1 c) 


completed by the perfect gas state relation. 


The geometry is made dimensionless by 
scaling the axial coordinate with the passage 
length L, the circumferential coordinate with the 
passage sector angle 4>, and the radial 
coordinate with the tip radius Rj. Let (R,x,0*) 
denote a dimensionless cylindrical coordinate 
system attached to the rotor. The x-axis lies 
along the axis of rotation and the rotating channel 
is restricted by the five walls: 0* = 0, 0* = 1 , R = 
Rh/Rj, R = 1 , and X = 1 . The equations are made 
dimensionless by introducing a reference length 
for the gradient operator, a reference time x, 
a reference density p„. and a reference speed of 
sound a_. As there is as yet no obvious way to 
assign definite values to the reference quantities, 
they will be left unspecified for now. 

Combinations of the reference scales are 
selected so that the following quantities, defined 
without asterisks, are dimensionless: 


P 


where h^ is the rothalpy, 
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and E is the total energy per unit mass in the 
rotating frame, 
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In addition to the ratios of geometric 
length scales (i.e., <I>Rt : R^ : L : LJ implicitly 
contained in the gradient operator, the equations 
Of nxDtion depend on three other dimensionless 
parameters 


B = 

T 
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The flow is modelled by the unsteady three- 
dimensional Euler equations [6] written in a 
rotating reference frame. 
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A proper scaling of the problem is required in 
order that the relative magnitude of each term in 


the equations of motion is indicated by a 
dimensionless factor preceding that term. Note 
that the wheel Mach number already indicates 
the n^gnitude of the centrifugal force but B and 
Rq demand further insights into the character of 
the problem. Boundary and initial conditions of 
the flow may be interrogated to define intrinsic 
reference scales of the problem. 

Boundary Conditions 

Boundary conditions at x=0 are 
formulated by introducing a gate function G(0) 
simulating the inlet-port opening and closing. An 
expression for the gate function is given as, 

G(e) = (ig) 


where H(0) is Heavyside’s unit step function, the 
angle 0 is an absolute angle measured in the 
frame of the stationary port, and 0p is the angular 
position of the port closing. A transformation from 
the absolute frame (0) to the rotating frame (0*) 
can be written as. 


0 = 0' + — f - —k Os0's1 , 
To ® 


where = 4>/i2 is the passage opening time. The 
above expression assumes that the origin of the 
two reference frames coincide at t=0. Introducing 
this transformation into equation (1g) gives an 
indication of whether or not a specific spatial 
location on the rotor passage at x=0 corresponds 
to a solid endwall { G (0*) = 0 ) or to the inflow- 
plane ( G (0*) = 1 ) at a given time. Inflow 
boundary conditions are discussed in the 
computational solution section. Free-slip wall 
velocity boundary conditions are: 


= 0 for. 


jf = 0 , X = 1 
, 0' = 0 , 0' = 1 


(1h) 


Initial Conditions 

At t=0, the passage is assumed to contain a gas 
at rest ( = 0 ) relative to the closed passage 

walls. Uniform thermodynamic conditions are 
prescribed (i.e., p = 1/y, p = 1). Although a radial 
variation of the thermodynamic properties can be 
prescribed by assuming homentropic solid-body 
rotation, this is probably not necessary for the 
range of parameters to be considered. 
Nevertheless, the present initialization procedure 
is thought to be close to the conditions in an 
actual rotor where the gas has been brought to 
rest through wave actions. 

Scalinq 

For the sake of sirrplicity, the uniform 
thermodynamic reference state is taken to be 
that which exists in the passage at the initial time 
(t = 0). In determining the reference length scale 
for the gradient operator and the reference time 
scale, attention should be focused on the region 
where and over what time period the dependent 
variables are changing most rapidly. It is 
anticipated that during the gradual opening 
period, most of the changes in the flow field will 
occur over a distance on the order of the 
passage width. This implies that 




= 




* Mq , f i 1 


Thus, during the gradual opening period, the 
wheel Mach number Mq is the time constant of 
the unsteady effects. That is, the higher the 
wheel Mach number, the more rapid will be the 
adjustment of the flow field in the passage to the 
application of the inflow conditions. 

Subsequent to the gradual opening 
period, it is anticipated that the length of the 
passage is the intrinsic length scale over which 
rapid changes will occur. Employing the 
longitudinal transport time Xl = l-/a,. as the 
intrinsic time scale results in the following 
dimensionless parameters 
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fl = = 1 . = (Z.//?;)Afo 

■t/. 



Hence, subsequent to the gradual opening 
period, unsteady and convective effects are of 
the same order of magnitude. 

It should be emphasized that the choice 
of intrinsic length and time scales depends on the 
parameter range under consideration. Observe 
that during the gradual opening period (i.e., 
physical times less than or equal to the opening 
time), a scaling is introduced which establishes 
as the dimensionless parameter indicating the 
relative order of magnitude of unsteady effects. 
However, subsequent to the gradual opening 
period, a different scaling, correctly indicating that 
unsteady and convective effects are of the same 
order of magnitude, should be adopted. 

For fixed geometry and boundary 
conditions, the governing equations are 
characterized by only one parameter, the wheel 
Mach number Mjj. The wheel Mach number has 
three different interpretations. First, it gives a 
measure of the relative strength of the centripetal 
acceleration. Second, during the gradual opening 
transient, the wheel Mach number can also be 
interpreted in terms of the ratio of the wave 
transport time across the width of the passage 
(4>RT/a_) to the passage opening time. Hence, for 
a fixed geometry and inlet conditions, increasing 
the wheel Mach number reduces the opening 
time relative to the wave transport time across 
the passage width. This implies a smaller 
distortion (i.e., relative to a planar orientation) of 
the wave fronts and also a reduction of the 
number of wave reflections. Third, M^ gives a 
measure of the deformation of the interface in 

terms of the ratio of tangential ( ~t ) to 
longitudinal ( ) extent of the interface 

during the gradual opening period. Moreover, the 
effects of Coriolis acceleration on the fluid contact 
interface will tend to be more pronounced due to 
a larger with increasing wheel Mach number. 


Computational Scheme 

The governing equations are solved 
using a four stage Runge-Kutta algorithm 
developed by Chima and Yokota [7]. This 
scheme employs explicit time-marching along 
with central-difference spatial discretization on a 
boundary-fitted H-mesh. The local truncation 
error is formally second-order both in space and 
time. Second and fourth difference artificial 
dissipation terms are used to capture 
discontinuities and to control non-linear 
instabilities. Further details concerning 
implementation of this scheme can be found in 
Reference 7. The treatment of the inflow 
boundary follows. 

The base flow at the inflow-plane (x=0) 
can be expressed as the sum of two 
components: a steady flow and an unsteady 
perturbatbn. This flow is assumed to evolve from 
the steady inlet-port flow, which is partially 
specified and partially computed from the solution 
interior to the rotor passage. It is this interior 
solution which is responsible for the unsteady 
perturbation. In addition, the base flow is 
assumed to be isentropic. The procedure 
adopted for describing the base flow at the 
inflow-plane consists of prescribing the incoming 
waves and using characteristic equations to 
supply the outgoing waves (see for example, 
Thompson [8] for more details). 


At each point on the inflow boundary four 
pieces of information must be furnished in order 
to calculate the flow variables. Assuming that the 
upstream axial Mach number is always subsonic, 
the Riemann invariant along the downstream 


running characteristic ( J* 




2a 

Y-1 


). the 


relative tangential velocity, the radial velocity, and 


the upstream entropy are specified. The Riemann 
invariant along the upstream running 

Oa 

characteristic ( J' = ) is extrapolated 

from the interior. This procedure allows all of the 
flow variables to be updated along the inflow 
boundary. 
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TABLE I.- SIMULATION PARAMETERS 


M„, M„ 


Rq ( based on Rt4> ) 

Rn (based on L) 

0.1 

0.573 

0.009 

0.150 

0.2 

0.287 

0.017 

0.300 

0.3 

0.191 

0.026 

0.450 

0.4 

0.144 

0.034 

0.600 

0.6 

0.096 

0.052 

0.900 


RESULTS AND DISCUSSION 

Five cases are presented to illustrate the 
flow dynamics of the gradual opening process 
and also the subsequent development of the 
interface during the time period for a single 
passage to traverse across an inlet port having 
an angular extent of approximately 13 passage 
widths (i.e., 6p = 14.2). A fixed passage geometry 
with a length-to-width ratio L/(Rt 4>) = 17.4, a hub- 
to-tip radius ratio R^/Rt = 0.934, and a length-to- 
tip radius ratio UR^ = 1.5 was used. This 
geometry reflects the NASA wave rotor 
experimental rig [9]. The uniform computational 
mesh consisted of 300 longitudinal nodes, 30 
radial nodes, and 30 circumferential nodes. A 
CFL = 0.5 was used to adequately resolve the 
gradual opening transient. 

Uniform inflow-plane conditions consisting 
of a downstream running Riemann invariant J* = 
7.08, upstream entropy P(/T/‘^’> = 0.126, and 
zero circumferential and radial relative velocity 
components (w^ = 0, Wr = 0) were assumed. 
These conditions are typical for a wave rotor 
passage filled with a quiescent charge of air at 
ambient conditions and encountering a perfectly 
matched (i.e., zero incidence to the rotor) inlet 
port supplying hot, high-pressure driver gas. The 
inflow-plane conditions are held fixed along with 
a specific heat ratio y = 1.4. The range of 
simulation parameters are given in Table I. Also 
listed in Table I are the ratios of full passage 
opening time (tJ to longitudinal transport time 
(\). 


Gradual Opening Transient Response 

The transient response of the flow field in 
the passage to the gradual opening process is 
illustrated in Figures 3, 4, and 5. These figures 
show pressure contours (Fig. 3), entropy 
contours (Fig. 4), and velocity vectors (Fig. 5) in 
a passage-to-passage plane (i.e., x-0^ at mid- 
height for several wheel Mach numbers and other 
parameters as given in Table I. Two instants 
during the gradual opening transient are 
presented, specifically the times when the 
passage is half opened (Fig. 3(a), 4(a), 5(a)) and 
fully opened (Fig. 3(b), 4(b), 5(b)). 

As the passage gradually opens, a series 
of curved pressure waves followed by a contact 
interface, as shown in Figures 3(a) and 4(a), 
emanate from the small opening. The waves 
penetrate into the passage and reflect from the 
trailing wall to interact with other incoming 
pressure waves and with the interface. Interaction 
of the reflected pressure waves with the interface 
produces baroclinic vorticity (i.e., x^p ) at 
the interface as is clearly evident in Figure 5. For 
Mo = 0.1 and 0.2, a large clockwise-oriented 
vortex is observed. A portion of the interface 
wraps around this vortex and assumes a spiral 
configuration. In this spiralling portion of the 
interface, the density gradient across the 
interface combines with the vortex-induced radial- 
pressure gradient (Fig. 3) to generate additional 
vorticity, possibly of opposite rotation. 

Note that the pressure waves undergo 
multiple reflections (Fig. 3) from both the trailing 
and leading walls of the passage. As the wheel 
Mach number is increased, a smaller number of 
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wave reflections occurs due to the more rapid 
adjustment of the flow. Thus, for increasing wheel 
Mach number during the gradual opening period, 
the baroclinic vorticity production at the interface 
due to interaction with reflected waves is reduced 
as is evident in Figure 5. 

It is observed from Figure 4 that the 
interface evolution during the gradual opening 
period exhibits increasing distortion. The interface 
distortion is caused by two mechanisms: (1) 
deformation of fluid elements comprising the 
interface as a result of gradual opening of the 
passage and (2) generation of vorticity (i.e., local 
rotation of fluid elements) by the actions of 
baroclinic torques, Coriolis forces, and vortex 
stretching. If the opening time were 
instantaneous, the interface would be a planar 
surface normal to the leading and trailing walls of 
the passage, but due to the finite opening-time, 
the interface deforms. 

The higher the wheel Mach number, the 
smaller the deformation as previously discussed. 
For Mo = 0.1 , a portion of the interface acquires 
a quasi-planar orientation (Fig. 4) extending 
across the passage width similar to that with 
instantaneous opening, followed by a trailing 
portion having a spiral configuration. The 
orientation of this spiral is consistent with the sign 
of baroclinic vorticity production (Fig. 5) due to 
wave reflections from the trailing wall. A similar 
but less developed interfacial structure is 
observed for Mo = 0.2. In comparison, the Mo = 
0.3-0.4 cases experience a much smaller 
deformation with no evidence of a spiral 
configuration. It can be inferred from the above 
considerations that whether the maximum 
transient response (i.e., the creation of a quasi- 
planar interface across the entire passage width) 
and the generation of a spiral configuration occur 
during or after the gradual opening period 
depends on the wheel Mach number. 

Since Rq (based on the passage width) 
is small (Table I) for the cases considered, the 
effect of Coriolis forces is negligible during the 
gradual opening period. Hence, the flow is mostly 
two-dimensional during this time period which 
implies that vortex stretching and tilting is small. 


Two-Dimensional Temporal Evolution 

The evolution of the interface and the 
propagation of compression or shock waves are 
shown in Figures 6(a)-(e) through density 
contours in a passage-to-passage plane at mid- 
height for wheel Mach numbers Mo = 0.1 , 0.2, 
0.3, 0.4, 0.6. Note that the dimensionless times 
are based on the longitudinal transport time scale 
(\) and the reference length is taken as the 
passage length L. Figures 6(a)-(e) spans the time 
period required for the passage to traverse the 
width of the inlet port. 

The plots suggest an interpretation of the 
interface in terms of two distinct structures: (1) a 
quasi-planar interface and (2) a trailing spiral 
interface resembling an evolving vortex sheet. 
The characteristics of the interlace evolution 
depend on two parameters: the wheel Mach 
number and Rq based on the passage length. 
For Mo = 0.1 -0.2 (Figs. 6(a)-(b)), the spiraling 
interface undergoes substantial distortion and 
eventually degenerates into a complex 
configuration (t = 1. 0-2.0). The quasi-planar 
interface exhibits a much smaller distortion with 
a slight bulge indicating a possible break-away 
from the spiral interface. 

For Mo = 0.4-0.6 (Figs. 6(d)-(e)), the 
spiral interface is smaller and better defined. The 
quasi-planar interface is more oblique and at t = 
1 .33, the interface appears to be spreading along 
both the trailing and leading walls. This is due to 
Coriolis and centripetal accelerations as will be 
shown later. For Mo = 0.3 (Fig. 6(c)) the interface 
orientation is intermediary to the Mo =0.1 and Mo 
= 0.6 cases. 

Also of interest in Figs. 6(a)-(e) is the 
formation of the shock wave and its interaction 
with the interface. It is apparent that at low wheel 
Mach numbers or large opening times, a 
relatively long time period is required for the 
formation of the shock wave. Figure 6(a) shows 
a series of compression waves propagating along 
the entire length of the passage for Mo = 0.1. 
Some of these waves reflect from the back wall 
and subsequently coalesce with part of the 
remaining incident compression waves to form a 
weak shock (t = 1.33). At a wheel Mach number 
of Mo = 0.6 (Fig. 6(e)), the shock is formed in a 
relatively short time period (t = 0.333). 
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The interaction of the shock with the 
interface is exhibited between t = 1 .667 and t = 
2.0. For Mo = 0.1, this interaction generates two 
distinct mixing pockets and at the same time, 
stretches the quasi-planar interface. Note also 
that the interaction of the shock with the interface 
spawns a transmitted shock along with a series 
of reflected pressure waves. An interesting result 
of the interface/reflected shock interaction is the 
counter-clockwise rotation of the quasi-planar 
interface for Mo = 0.1 -0.3, whereas for Mo = 0.4- 
0.6, a complex interfacial pattern develops. 

Further insight into the evolution of the 
interface is given in Figure 7 which shows 
contours of the radial component of vorticrty in a 
passage-to-passage plane at mid-height for Mo 
= 0.1, 0.3, 0.6 at t = 0.667, 1.33, and 2.0. A 
negative vorticity implies a clockwise rotation. 
The vorticity distributions are consistent with the 
interface orientations in Figs. 6 (a)-(e). For Mo = 
0.1 at t = 0.667 (Fig. 7(a)), the negative vorticity, 
generated by wave reflections from the trailing 
wall, appears to be surrounded by a region of 
positive (counter-clockwise) vorticity, possibly 
generated by wave reflections from the leading 
wall. There is a local positive leading-wall vortex 
for Mo = 0.6 at t = 0.667. Also, at t = 1 .33 for Mo 
= 0.6, the negative spiral vortex is sandwiched 
between positive vortices at both walls. 

Interaction of the reflected shock with the 
interface results in a complex vorticity distribution 
as is evident in Fig. 7 at t = 2.0. Multiple pockets 
of negative and positive radial vorticity are 
generated for Mo = 0.1. For Mo = 0.6, the 
negative spiral vorticity is not discernible. 

Three-Dimensional Effects of Passage Rotation 

The three-dimensionality of the interface 
development is illustrated in Figure 8 at two 
instants subsequent to the gradual opening 
transient. Entropy contours on passage-to- 
passage planes at hub, mean, and tip radii are 
shown for wheel Mach numbers Mo = 0.1, 0.3, 
0 . 6 . 

For Mo = 0.1, there does not appear to 
be any substantial three-dimensionality in the 
interlace configuration. Similar patterns can be 
found at hub, mean, and tip. However, for Mo = 
0.3, there is more intense mixing at the tip spiral 


vortex relative to the hub. The center of the spiral 
vortex is at roughly the same axial position for 
hub, mean, and tip sections. For Mo = 0.6, a 
radial distortion of the interface with driven fluid 
(lower entropy) on top and driver fluid (higher 
entropy) on the bottom is observed. This is a 
result of baroclinic vorticity generation by the 
centripetal acceleration at the interface. 

For Mo = 0.6, there exists circulatory 
flows normal to the mean flow due to vorticity 
generation by Coriolis forces. This is 
demonstrated by the presence of interfacial 
material along the trailing and leading walls of the 
passage at mid-height for Mo = 0.6. Material from 
hub and tip migrate spanwise towards mid- 
height. 

Wave Rotor Design Implications 

The static pressure and axial velocity 
distribution along the passage length at mid- 
height and mid-passage are shown in Figure 9 
for t = 0.667. This figure gives a global 
perspective of the effects of various gradual 
opening times. Also included in Figure 9 is the 
instantaneous opening case as a reference. It 
can be seen that reducing the gradual opening 
time (i.e., irrcreasing MJ for a fixed geometry 
allows a more rapid shock formation and 
propagation. The time for shock formation is 
significantly greater than that with instantaneous 
passage opening. This has strong implications for 
port timing and wave arrangement. 

Increasing the opening time (i.e., 
reducing MJ generates substantial blockage near 
the interface resulting in a locally higher 
throughflow velocity at a reduced static pressure. 
This is a loss mechanism (i.e., mixing loss 
between driver and driven gas) that should be 
addressed by wave rotor designers. The flow 
blockage is caused by the vortical structures that 
are generated by the gradual opening process 
and sustained by the action of centripetal and 
Coriolis accelerations at the interface. For a fixed 
overall geometry, it may be possible to reduce 
this blockage by profiling the passage walls near 
the inlet. This profiling needs to be performed in 
a three-dimensional fashion if the adverse effects 
of the centripetal and Coriolis accelerations are to 
be countered. It should also be noted that an 
optimum wheel Mach number allowing a 
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compromise between the two-dimensional 
adverse effects of gradual opening and the 
associated three-dimensional effects resulting 
from Coriolis and centripetal accelerations may 
be found. 


CONCLUSIONS 

A numerical simulation of the flow 
dynamics of gradual opening of a wave rotor 
passage shows that the flow development is 
complex and three-dimensional. For fixed 
geometry and boundary conditions, the response 
of the flow in the rotor passages depends on the 
wheel Mach number. At low wheel Mach 
numbers (i.e., Mo < 0.3) or long opening times, 
the interface between in-flowing driver gas and 
the driven gas in the passage becomes very 
distorted due to interactions with multiple 
reflected compression waves from the trailing and 
leading walls of the passage. Vortices are 
produced at the interface and it is observed that 
these vortices can serve as a stirring mechanism 
to promote mixing between driver and driven fluid 
along the interface. For higher wheel Mach 
numbers (i.e.. Mo > 0.2) or short opening times, 
the amount of compression wave reflection is 
reduced which results in a smaller distortion of 
the interface due to interactions with reflected 
compression waves. However, the higher 
centripetal and Coriolis accelerations generate 
additional vorticity which causes a three- 
dimensional interface distortion. 

The development of the wave system is 
also dependent on the wheel Mach number or 
passage opening time. For low wheel Mach 
numbers or longer opening times, a substantial 
time period is required for the formation of 
shocks. Due to the highly distorted interface, 
subsequent interaction of the reflected shock with 
the interface results in the creation of a complex 
flow pattern. 

It can be concluded that the gradual 
opening process and the influence of rotation on 
the flow patterns are intimately coupled. This 
should be kept in mind during the wave rotor 
design process, not only for proper port-timing 
and wave arrangement but also because of the 
losses which will occur due to mixing and wave 


reflections. Rotating faster in order to reduce the 
adverse effects of gradual opening may create 
other adverse effects due to Coriolis and 
centripetal accelerations. The results presented 
herein and additional results with the inclusion of 
viscous forces may help guide the designer in 
achieving an optimum compromise. 
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Figure 1. Wave Rotor Schematic Diagram. 



Figure 2. Geometric configuration for simulating the gradual opening process in a wave rotor. 
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o L 



(a) half opened t = 0.5 



Figure 3. Pressure contours in a passage-to-passage plane at mid-height illustrating the evolution of compression 

waves during the gradual opening transient for various wheel Mach numbers (Mo) and ratios of opening time to 

longitudinal transport time (t /t. ). 
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Figure 4. Evolution of the interface for different wheel Mach numbers (Mo) at two instants during the 
gradual opening transient. The contours correspond to entropy In a passage-to-passage plane at mid-helght. 
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(a) half opened t = 0.5 






(b) fully opened t = 1.0 

Figure 5. Velocity vectors in a passage-to-passage plane at mid-height during the gradual opening transient. 
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t=.287 t=.333 t=.573 t=.667 t=1.00 



Figure 6(a). Density contours in a passage-to-passage plane at mid-height for a wheel Mach number Mo=0.1. 
CW-compression wave; RCW-reflected compression wave; RP-reflected pressure wave; TS transmitted shock. 
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Figure 6(b) Density contours in a passage-to-passage plane at mid-height for a wheel Mach number Mo=0.2. 
CW-compression wave; RS-reflected shock: RP-reflected pressure wave. 
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Fiigure 6(c). Density contours in a passage-to-passage plane at mid-height for a wheel Mach number Mo=0.3. 
S-shock: RS-reflected shock: RP-reflected pressure wave. 


21 


Page intentionally left blank 



t=.072 t=.144 t=.333 


t=.667 


t=1.00 t=1.333 t=1.667 


t=2.0 



Figure 6(d). Density contours in a passage-to-passage plane at mid-height for a wheel Mach number Mo=0.4. 
S-shock; RS-reflected shock: RP-reflected pressure wave. 
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Figure 6(e). Density contours in a passage-to-passage plane at mid-height for a wheel Mach number Mo=0.6. 
S-shock; RS-reflected shock; RP-reflected pressure wave; TS-transmitted shock. 
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(a) t = 0.667 




r 

(c) t = 2.0 

Figure 7. Contours of the radial component of vorticity in a passage-to-passage plane at mid-height. 
(-) corresponds to a clockwise rotation; (+) corresponds to a counter-clockwise rotation. 
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Mo = 0.1 



Mo = 0.3 ^ driven fluid (low entropy) 



X 


(a) t = 1.0 


spiral , vortex 



driver fluid (high entropy) 


(b) t = 1.33 

Figure 8. Entropy contours for various wheel Mach numbers at two instants subsequent to the gradual 
opening transient. The three-dimensionality of the interfacial structure is illustrated for Mo = 0.3 and 0.6. 
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X 

(a) static pressure 



(b) axial velocity 

Figure 9. Static pressure and axial velocity distribution along passage length at mid-height and mid-passage 
for t = 0.667. 
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